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Abstract 

The chemical and photometric evolution of star form- 
di&k salaKies is investigated, Numerical simulatic 



mergers (Mihos & Hernquist 1996) and of galaxy forma- 
tion and evolution (Katz 1992). A second good example is 
an GRAPE - SPH c ode ( |Steinmetz fc MuUer 199^ ; ^tchT 



r .-i 1 1 • ifl u J metz & Muller 1995) which was successfully used to model 
01 the complex gasdjrnamicai nows are based on our own . — f 



coding of the Chemo - Dynamical Smoothed Particle Hy- 
drodynamical (CD - SPH) approach, which incorporates 
the effects of star formation. As a first application, the 
model is used to describe the chemical and photometric 
evolution of a disk galaxy like the Milky Way. 

Key words: star formation: chemical and photometric evo- 
lution: SPH code 



1. Introduction 

Galaxy formation is a highly complex subject requiring 
many different approaches of investigation. Recent ad- 
vances in computer technology and numerical methods 
have allowed detailed modeling of baryonic matter dy- 
namics in a universe dominated by coUisionless dark mat- 
ter and, therefore, the detailed gravitational and hydro- 
dynamical description of galaxy formation and evolution. 
The most sophisticated models include radiative processes. 



the evolution of disk galaxy structure and kinematics. 



star formation and supernova feedback ( e.g. Katz 1992; 
Steinmetz &: Muller 1994|; [Friedh fc Benz 1995). 



The results of numerical simulations are fundamentally 
affected by the star formation algorithm incorporated into 
modeling techniques. Yet star formation and related pro- 
cesses are still not well understood on either small or large 
spatial scales. Therefore the star formation algorithm by 
which gas is converted into stars can only be based on sim- 
ple theoretical assumptions or on empirical observations 
of nearby galaxies. 

Among the numerous methods developed for modeling 
complex three dimensional hydrodynamical phenomena. 
Smoothed Particle Hydrodynamics (SPH) is one of the 
most popular ( [Monaghan 1992 ). Its Lagrangian nature 
allows easy combination with fast N - body algorithms, 
making possible the simultaneous description of complex 



gas-stellar dynamical systems (Friedli & Benz 1995). As 



an example of such a combination, an TREE - S PH code 
([Hernquist fc Katz 19891 ; [Navarro fc White 199 j ) was suc- 
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Figure 1. Luminosity evolution of the model galaxy. 



2. Model 

The hydrodynamical simulations are based on our own 
coding of the Chemo - Dynamical Smoothed Particle Hy- 
drodynamics (CD - SPH) approach, including feedback 
through star formation (SF). The dynamics of the "star" 
component is treated in the frame of a standard N - body 
approach. Thus, the galaxy consists of "gas" and "star" 
particles. For a detailed description of the CD - SPH code 
(the star formation algorithm, the SNII, SNIa and PN 
production, the chemical enrichment and the initial condi- 



tions) the reade r is referred to Berczik fc Kravchuk (1996) 



cessfully applied to the detailed modeling of disk galaxy of our algorithm. 



Berczik (1999) . Here we briefly decsribe the basic features 



Proc. of 33rd ESLAB Symp. "Star formation from the small to the large scale" (F. Favata, A. A. Kaas & A. Wilson eds., ESA 
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The chosen "gas" particle produces stars only if the 
above condition holds over the time interval exceeding its 
free - fall time: 



(2) 

'gas" particles remain cool, 



<tff- We rewrite these conditions following Navarro 
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(3) 

We set the value of pcrit — 0.03 cm~^. 
When the collapsing particle i is defined, we create the 
new "star" particle with mass m'**"'' and update the "gas" 
particle rrii using these simple equations: 



0.0 3.0 6.0 9.0 12.0 15.0 

t (Gyr) 

Figure 2. Photometric evolution of the model galaxy. 



mj = (1 - e) • mj. 



(4) 
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In the Galaxy, on the scale of giant molecular clouds, 
the typical values for SF efficiency are in the range e ~ 
0.01 ^ 0.4 (puerr et al. 1982|; [Wilking fc Lada 1983|). 
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I '■""■■-•"■'Wff'iii, 



We did not fixe this value but rather also derived e 
from the "energetics" condition: 



e = 1 - 



th 



E: 



gr 



(5) 



At the moment of birth, the positions and velocities of 
new "star" particles are set equal to those of parent "gas" 
particles. Thereafter these "star" particles interact with 
other "gas" and "star" or "dark matter" particles only by 
gravity. 

For the thermal b udget of the ISM, SNIIs play the 
main role. Following to |Katz (1992] ; |Friedh fc Benz (1995) 



we assume that the energy from the explosion is converted 
totally to thermal energy. The total energy released by 
SNII explosions (10'''' J per SNII) within "star" parti- 
cles is calculated at each time step and distributed uni- 



formly between the surrounding "gas" particles (Raiteri 
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Figure 3. Color index evolution of the model galaxy. 



We modify the standard SPH SF algorithm (|Katz 1992| ; 
Stcinmctz fc MuUcr 199^ ; [Steinmetz fc MuUer 1995D , tak- 
ing into account the presence of chaotic motion in the 
gaseous environment and the time lag between the initial 
development of suitable conditions for SF, and SF itself. 

Inside a "gas" particle, the SF can start if the absolute 
value of the " gas" particle gravitational energy exceeds the 
sum of its thermal energy and energy of chaotic motions: 



Ef ||> -f K 



In our SF scheme, every new "star" particle represents 
a separate, gravitationally bound, star formation macro 
region (like a globular clusters). The "star" particle has its 
own time of birth tbegs f which is set equal to the moment 
the particle is formed. After the formation, these particles 
return the chemically enriched gas into surrounding " gas" 
particles due to SNII, SNIa and PN events. 

We concentrate our treatment only on the production 
of "'^^O and ^^Fe, yet attempt to describe the full galactic 
time evolution of these elements, from the beginning up 
to present time ( i.e. t^voi ~ 15.0 Gyr). 

The code also includes the photometric evolution of 
each "star" particle, b ased on the idea of the Single Stel- 
lar Popu lation (SSP) ( Bressan et al. 1994 ; Tantalo et 
1996| ). 
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Figure 4- Evolution of the model galaxy in the (B-V) vs. V Figured. Evolution of the model galaxy in the (U-B) vs. 



plane. 



(V-K) plane. 
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Figured. Evolution of the model galaxy in the (U-B) vs. 
(B- V) plane. 



At each time - step, absolute magnitudes: M(7, Mb, 
My, Mr, M/, Mi<-, Ma/ and MboZ are defined separately 
for each "star" particle. The SSP integrated colours (UB- 
VRIKM) are taken from [Tantalo et al. (1996)} The 
spectro - photometric evolution of the overall ensemble 
of " star" particles forms the Spectral Energy Distribution 
(SED) of the galaxy. 



We do not model the energy distribution in spectral 
lines nor the scattered light by dust. However according 



to Tantalo et al. (1996) our approximation is reasonable, 



especially in the UBV spectral brand. 



3. Results 

The model presented descibes well the time evolution of 
the basic chemical and photometric parameters of a disk 
galaxy similar to the Milky Way. The metallicity, luminos- 
ity and colors obtained are typical of such disk galaxies. 

— Figure 1 . Luminosity evolution of the model galaxy. 

— Figure 2 . Photometric evolution of the model galaxy. 

— Figure 3. Color index evolution of the model galaxy. 

— Figure i. Evolution of the model galaxy in the (B-V) 
vs. V plane. 

— Figure |[ Evolution of the model galaxy in the (U-B) 
vs. (B-V) plane. 

— Figure ^ Evolution of the model galaxy in the (U-B) 
vs. (V-K) plane. 

— Figure 0. Evolution of the model galaxy in the [Fe/H] 
vs. V plane. 

— Figure^. The "real" Mgtar/Ly evolution of the model 
galaxy 
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Figure 7. Evolution of the model galaxy in the [Fe/Hj vs. 
V plane. 



1.0 




t (Gyr) 

Figures. The "real" Mstar/Lv evolution of the model 
galaxy. 
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